Shear and bulk viscosities for pure glue matter 
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Shear 77 and bulk ( viscosities are calculated in a quasiparticle model within a relaxation-time ap- 
proximation for pure gluon matter. Below T c , the confined sector is described within a quasiparticle 
glueball model. The constructed equation of state reproduces the first-order phase transition for the 
glue matter. It is shown that with this equation of state, it is possible to describe the temperature 
dependence of the shear viscosity to entropy ratio 77/s and the bulk viscosity to entropy ratio f/s in 
reasonable agreement with available lattice data, but absolute values of the £/s ratio underestimate 
the upper limits of this ratio in the lattice measurements typically by an order of magnitude. 
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I. INTRODUCTION 

The high-energy heavy-ion collisions at Super Proton 
Synchrotron (SPS) and Relativistic Heavy-Ion Collider 
(RHIC) energies have shown evidence of a new state of 
matter characterized by very low shear viscosity to en- 
tropy density ratio 77/s similar to a nearly ideal fluid 0- 
4]- Lattice calculations indicate that the crossover re- 
gion between hadron and quark-gluon matter has been 
reached in these experiments. On the other hand, lat- 
tice calculations performed in gluodynamics (GD) clearly 
demonstrate the occurrence of there occurs the first-order 
phase transition. 

The shear 77 and bulk £ viscosities are parameters that 
quantify dissipative processes in the hydrodynamic evolu- 
tion of a fluid. It is known that the behavior of transport 
coefficients is sensitive to the presence of phase transi- 
tions in a medium (see 0-[LH and references therein). 
Values of the bulk and shear viscosities near the phase- 
transition critical temperature T c affect the hydrody- 
namic evolution of the medium and may influence ob- 
servables. 

Lattice quantum chromodynamics (QCD) is the most 
powerful technique to extract nonpcrturbative informa- 
tion on an equation of state (EoS). When experimental 
data are lacking, lattice data are often used to fit model 
parameters. For pure gluon SU(3) theory, the EoS was 
computed on the lattice more than a decade ago [l^l ■ Re- 
cently, much more accurate data have been obtained fl3j . 

Among various existing phcnomenological approaches, 
quasiparticle (QP) models are used to reproduce results 
obtained in the lattice QCD. In the case of GD, the QP 
models rely on the assumption that for the temperature 
T above the critical one, T > T c , the system consists 
of a gas of massive deconfined gluons. In the confined 
phase, at T < T c , the glue matter is considered as a gas 
of massive glueballs. 

In this paper we aim to investigate the behavior of 
viscosity coefficients for a gluon system that exhibits a 
deconfincment phase transition. The phenomenological 



QP model is applied to describe available lattice data on 
the EoS. Shear and bulk viscosities are calculated within 
a relaxation-time approximation. 



II. EQUATION OF STATE OF GLUE MATTER 

In the QP approach the system of interacting gluons 
is treated as a gas of noninteracting quasiparticlcs with 
an effective mass m g (T), which depends on T as Q 



m 2 g (T) = ^ g\T) T 



(1) 



with the temperature-dependent strong interaction con- 
stant 
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where parameters T s /T c =0.5853, A =3.3 are taken to fit 
the new lattice data, see below, and a number of colors 
N c =3. The energy density and the pressure acquire then 
the following forms: 



= el d (T,m g (T)) + B(T), 



P g (T) = & f °° P 2 dp p 4 mp(E ' /TH - B(T) 
= Pl d {T,m g {T))-B(T), 



(3) 



(4) 



where the degeneracy factor d g = 2(N^ — I) = 16 for the 
SU(3) gluodynamics, e ld and P g d are the energy density 
and the pressure of the ideal gas of massive gluons. The 
temperature-dependent function B{T) in Eg. (j5t results 
from the thermodynamical identity, sec Ref. |14| . 



T^-P(T) = e(T), 



(5) 
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which leads to the equation for B(T): 



dB{T) 
dT 



ef - 3Pf dm g 



m„ 



dT 



(6) 



Dealing only with gluon degrees of freedom, one as- 
sumes that the matter at T < T c (the "hadronic" phase) 
consists of glueballs. While the meson scattering ampli- 
tude is parametrically suppressed as 1/N C , the scatter- 
ing amplitude between glueballs scales as 1/N% [l5[ and 
therefore the system can be considered as a nonintcract- 
ing Bose gas of glueballs. Expected glueball masses are 
high, rrigb >1 GeV, and thereby only the lowest- lying 
glueball states contribute to the EoS at the tempera- 
tures of our interest. It is difficult to single out which 
states of the observed hadronic spectrum are glueballs 
because of a lack of knowledge of decay properties and 
the existence of a strong mixing between glueballs and 
quark states (l6| . However, using typical constant values 
for the lowest-lying glueball masses within a statistical 
model, one fails to reproduce the strong increase of ther- 
modynamical variables near T c [TtJ ■ The T-behavior of 
masses for two lowest-lying scalar ++ and tensor 2 ++ 
glueballs was investigated on the lattice in [TJJ. There- 
fore, below we follow the SU(3) lattice GD results. It 
was shown that the pole mass m g b(T), the Breit-Wigner 
mass rh g b(T), and the thermal width T g b are linked as 
follows : 



m gb (T) » m gb (T) - 2T + y£p - T* gb (T) 



(7) 



With the help of the Ansatz m g b(T) = m® b , i.e., that the 
glueball Breit-Wigner masses are given by the particle 
data group (PDG) values, 

r g b = b gb (T-T gb ) e(T-T gb ) for T gb < T < T c (8) 

and recommended parameters b g b(0 ++ ) =4.23 and 
b g b(2 ++ ) =7.152, the relation (J7]) reproduces quite well 
the lattice results in the measured range 0.5T C < T < 
T c = 265 MeV 0]. In our consideration, we limit our- 
selves to the two above-mentioned species of glueballs, 
the only ones for which lattice data are available. 

With the temperature-dependent glueball masses, a 
statistical treatment of glueballs needs an additional re- 
quirement of thermodynamic consistency. It has been 
satisfied in the same way as outlined above for gluons; 
see Eqs. ([5]) and (|6]) above. 

To describe glue matter in the whole range of temper- 
atures, we use the first-order phase-transition model in 
accordance with lattice results for GD. Thus, one should 
conjugate the pure gluon (g) and the glueball (gb) phases 
by making use of the Gibbs conditions at the transition: 



Tf. = Tf = T r 



P g {T c ) = P gb (T c ) 



(9) 



We use the value T c =265 MeV for the first-order phase 
transition, in agreement with the lattice SU(3) GD [l2|, 
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FIG. 1: The reduced pressure (a) and the energy density (b) 
of the glue matter. The solid line corresponds to taking into 
account two glueballs, scalar ++ with m° gb = 1470 MeV and 
tensor 2 ++ with m° gb = 2150 MeV. Experimental points are 
the old Karsch's (filled squares) and the new Panero's 
(circles) QJ] lattice results. The region near T c is enlarged in 
the inset. 




FIG. 2: The trace anomaly (a) and the reduced enthalpy or 
entropy (b) of the glue matter. All notations are the same as 
in Fig. [1] 



In Fig. [TJ we compare the model results for the pres- 
sure and the energy density with the lattice data. Values 
arc normalized to those in the Stcfan-Boltzmann (SB) 
limit. Both e/esB and P/Psb increase fast and mono- 
tonically with the temperature above T c , but, as we see, 
the Stefan-Boltzmann limit is not saturated up to 3T C . 
Two sets of lattice data are similar qualitatively, but old 
data [2(| are appreciably higher than the new ones [Jj| at 
T > 1.5T C . Note that the new lattice data are extended to 
the region of the glueball phase, T < T c . Our QP model 
is in a reasonable agreement with the new lattice data, 
except for the narrow vicinity to the left of T c , where 
model predictions are evidently below the lattice points 
[see insertions in Fig. [TJa) and Fig. Gib)]. Due to large 
glueball masses, this result is not changed if one adds the 
next two to three glueball states to our two lowest-lying 
glueball states, although there exist statements claiming 
that a good agreement with lattice data near T c can be 
reached only if the whole high-lying glueball spectrum 
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of the Hagedorn-type 21] or glueball condensate UJ is 
additionally included. 

The interaction measure or trace anomaly (e — 3P)/T 4 
and the reduced enthalpy (e + P)/(esb + -Psb) arc demon- 
strated in Fig.[5{a) and Fig.[3{b), respectively. The pres- 
ence of a nonzero trace for the energy-momentum tensor 
relates to the breaking of the scale and conformal in- 
variancc. Again, a nice overall agreement is observed 
between the QP model and the new set of the lattice 
data [Hj], including a region near T c , for T > T c . The 
reduced enthalpy in Fig. [2jb) for the pure gluon system 
is just the reduced entropy, s/ssb, which is thereby also 
reproduced by our QP model. Thus, we see that the 
developed QP model successfully describes the thermo- 
dynamic properties of the glue matter. 



III. CALCULATION OF VISCOSITY 
COEFFICIENTS 

In principle, it is possible to compute the shear and 
bulk viscosities directly from GD at finite temperature 
using Kubo formulas. However, in practice, this is quite 
difficult because GD is generally a strongly interacting 
theory with an unknown mechanism of the confinement. 
Essential assumptions of our kinetic approach are that 
quasiparticles are well defined, elementary interactions 
are local, and the dynamics can be described in the 
relaxation-time approximation. 

Derivation of viscosity coefficients starts with the ex- 
pression for the energy-momentum tensor for quasi-free 
[43f boson quasiparticles of spices a : 



(10) 



dT 



d 3 p a 
(2tt) 3 



{E a (p a ,r),p a ) 



where d a is the degeneracy factor. The QP distribution 
function F a fulfills the QP kinetic equation. Wc assume 
that gluon and glueball masses are given by Eqs. ([1]) and 
0, respectively. The QP energy is determined by 



E a (p) = y/p2 + ml(T,F a ) 



(11) 



Below, we consider only collisional sources of the vis- 
cosity. Applying the relaxation time approximation to 
the relativistic QP kinetic equation, we arrive at the ex- 
pression for the variation of the energy-momentum tensor 
(flOl) near the local equilibrium state: 



(12) 



loc.eq. 



where r a denotes the relaxation time of the given species, 
which generally depends on the QP momentum p a . The 



local equilibrium distribution function for a boson is as 
follows: 



F^(p a ,x a ) 



1 



(13) 



~ (l,u) for \u\ <C 1. Performing a variation in (|12l) . 
we did not vary quantities that may depend on the distri- 
bution function only implicitly, such as E a , since by doing 
this one may arrive at the relaxation-time form of the QP 
kinetic equation. Besides, in the gluon-glueball model 
used here, only equilibrium values m a (T) are known and 
we are actually not able to find SE a [F] . 

The shear and bulk viscosities can be expressed 
through the variation of the energy-momentum tensor 
as follows: 



STi, 



-C <%V -u — r) Wij, 



(14) 



with Wki = dkUi + diu k - -S k i dtu 1 



Here and below, Latin indices run 1,2,3. To find the 
shear viscosity, we set i ^ j in (|14[) . To find the 
bulk viscosity, we substitute i = j in (|14j) and use that 

■^loc.eq = 3-Ploc.cq- 

Taking derivatives dF^ c eq - / dx^ 1 in Eq. (fL2|) and using 
(fl"4| as a definition of viscosity coefficients, by straight- 
forward calculations we find expressions (see [23 - I25T ]) for 
the shear viscosity, 



1 



dTT a ^F^ (1 + FD 



(15) 



and for the bulk viscosity [4j 



1 



dTr a ^F^ {l + F?)Q a , (16) 



3T ^ 

a 

where the EOS-dependent Q a factor is given by 

Ma 



Pi 



E a -T- 



dT 



(17) 



and cj = is the speed of sound squared. 

Simplifying, instead of the momentum-dependent 
value r a , one may use the averaged partial relaxation 
time f a related to the cross section as 



(18) 



where n a > is the particle density of a species, er^, = 

/ dcosO da(aa —> aa )/dcosd (1 — cos 9) is the trans- 
port cross section, in general, accounting for in-medium 
effects, and v aa > is the relative velocity of two colliding 
particles a and a in the case of binary collisions. Angu- 
lar brackets denote a quantum-mechanical statistical av- 
erage over an equilibrated system. However, one should 
bear in mind that averaged values f" 1 given by Eq. (|18[) 
yield only a rough estimate for the values r" 1 . 
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IV. RESULTS FOR VISCOSITIES 

Below, the shear and bulk viscosities are calculated 
with the help of Eqs. (|15|) and ([To]) , respectively. The 
only quantity that should still be specified is the relax- 
ation time r a . 

Calculations of the relaxation time f a of partons al- 
ready in the lowest order in the running coupling constant 
g 2 require summation of infinitely many diagrams. Re- 
summation of the hard thermal loops results in the width 
f -1 of partons ~ g 2 Tln(l/g) [26[. Based on this fact, the 
following parametrization was used for gluons [H, H3] : 



= N r 



9 * T 



-in 



9' 



(19) 



with the strong interaction coupling constant (|2|) and a 
tunning parameter c. The relaxation time for a mixture 
of scalar and tensor glueballs was estimated according to 
Eq. (|18p assuming the glueball scattering cross section 
(Tgi, =30 mb to be isotropic. 
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FIG. 3: The ratio of the shear viscosity to entropy density 
for a pure glue matter. Solid and dashed lines are for our QP 
two-phase gluon-glueball model with two different choices of 
the coefficient c in Eq. (|19[) ; see the text for details. The 
vertical short-dashed line joins two boundary points of the 
mixed gluon-glueball phase. The dot-dashed line shows vis- 
cosity calculations with the relaxation time f 3 given by Eq. 
(|20l) . The horizontal dotted line is the rj/s = 1/4-7T bound. 
The lattice gauge SU(3) data with 16 3 ■ 8 and 24 3 • 8 lat- 
tice are from Refs. [2S| (triangles and squares) and [53] (filled 
circles). The shaded region corresponds to the perturbative 
result (cited from [2^1. 

Comparison between the GD lattice data [29|, H(| and 
our QP results for the shear viscosity of the glue matter 
is presented in Fig. [3] The magnitude of the rj/s ratio in 
our model is defined mainly by the value of the relaxation 
time (flU)) . The solid line shows results of our calculation 
provided we use the recommended value c =14.4, though 
our parameters of Eq. ([2]) are slightly different from those 
used in Ref . [J • This c- value was tuned in Re f. [3 to the 
old lattice data for thermodynamic quantities [12[ • As we 
see, the rj/s ratio gets discontinuity at T = T c with more 



than by an order of magnitude lower value at T — > T c + 
(in gluon phase) than at T — > T c — (in the glueball 
phase). Also, the solid curve lies reasonably close to the 
points (filled circles, triangles, and squares) and its value 
in the minimum is slightly below the AdS / CFT (Anti de 
Sitter/Conformal Field Theory) 1/47T bound [31( (com- 
pare with the dotted curve). Preserving the form of the 
relaxation time (|19p . we can still increase rj/s values by 
tuning the parameter c. Taking c = 11.44, wc achieve the 
limit case f g ~ l — > for T — > T c + (full transparency). 
In this case (see the dashed line in Fig. [3] for T > T c ), 
we may reach a sli ght ly better overall agreement with 
the lattice data [29L l30j. and the 1/4-7T bound is achieved 
at the minimum. Varying the c-value in the interval 
11.44 < c < 14.4, one may simulate different values of 
the rj/s jump at T = T c , but for temperatures T > 1.5T C 
the rj/s ratio changes only slightly, demonstrating a slow 
increase with the growing temperature. Thus, bearing in 
mind large error bars in the lattice data, we are able to 
conclude that the results of the two-phase gluon-glueball 
model developed here are consistent with the existing 
lattice results [13, [3(| . The perturbative regime (see the 
shaded region) is not achieved up to very high tempera- 
tures. 

The rj/s ratio for the pure gluon phase in the range of 
T ~ (1 — 2)T C was also evaluated in Ref. [32j|. The model 
employs the QP Ansatz for EOS successfully tested to 
describe old lattice results [28|. In Ref. [32j], viscosity is 
treated by means of the kinetic theory for gluon quasipar- 
ticles. It is of interest that the model, being consistent 
with the old (and less accurate) lattice data for viscos- 
ity (28| and thermodynamics [2(| that overestimate pres- 
sure at T > 1.5T C (as follows from the comparison with 
new data, see Fig. []}, predicts a stronger temperature 
dependence of rj/s at T > T c than our model, which in 
turn is consistent with the new lattice data [23, H3| ■ The 
crucial point here is that the gluon relaxation time is 
defined essentially differently: 



fr 1 



where 



BKR 



=6. 



a„/(327r 2 )T g 4 log( S 7r/ 5 2 ) 



(20) 



Here f BKR 



cx (? 4 , coming from the 



a s ~ g dependence of the gluon-gluon transport cross 



section cr* , in the relaxation-time expression (fl8|) , as was 
estimated in the early work of Hosoya and Kaiantie [3(| , 
whereas the above used t^ 1 oc g 2 |4|. |27|. 

The rj/s ratio obtained with the relaxation time ([2TJ| 
is plotted in Fig. [3] by the dot-dashed line. This re- 
sult for T < 2T C recovers that of Ref. [13], but it differs 
significantly from those calculated with Eq. (jT9")) . Us- 
ing recent lattice results for higher T (29[, let us try to 
make a choice between two paramctrizations of relaxation 
times (3U and ([2D]). We observe that for T > 10T C , shear 
viscosity calculations with ([2TJ]) demonstrate a noticeable 
growth exceeding lattice data and even a perturbative 
estimate (v/s) pcrt « 0.8 - 1.0 for T ~ (2 - 20)T C . The 
perturbative result is taken from Fig. 2 of [Hj] (77 was cal- 
culated in [33j] and s in [HI). In contrast, predictions of 
our QP model with relaxation time (|19[) arc in reasonable 
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agreement with the lattice results and do not contradict 
perturbative estimates. Recently, the relaxation time f g 
was estimated in Ref. (HJ according to Eq. (fT8|) from an 
analysis of cross sections of the gg — > gg and gg — > ggg 
processes. It was found that ?7/s=0.13 and 0.076 for val- 
ues a s =0.3 and 0.6, respectively [which correspond to 
temperatures T/T c =2.6 and 1.36 provided Eq. ([2j is 
used] . If these points were plotted in Fig. [31 they would 
be quite consistent with our QP model results. This can 
be considered as an additional numerical argument in fa- 
vor of using Eq. (|19p as a rather appropriate phenomeno- 
logical expression. However, we point out that the above 
arguments do not allow us to choose between functional 
dependencies (fT9]) and (|20|) . One could also use expres- 



sion pop with a smaller coefficient a v , choosing the latter 
value to fit the lattice results for the r\j ' s ratio. Since un- 
certainties in the modern lattice data are large, we will 
not do such an additional tuning restricted by the results 
that we have demonstrated in Fig. [3] 

The measured lattice points for the ratio of the bulk 
viscosity to the entropy density are plotted in Fig. 2] to- 
gether with different model results. In a broad range of 
temperatures, the global behavior of lattice data can be 
roughly approximated as C/s = 0.02/ ^T/T c — 1 (see the 
short-dashed curve). The reduced bulk viscosity (/s cal- 
culated in our two-phase gluon-glueball model following 
Eqs. ([TBI and (|19[) is shown by the solid line for c = 14.4 
and by the dashed line for f^ 1 vanishing at T c . Val- 
ues of C/s for both curves noticeably underestimate the 
corresponding values on the approximating short-dashed 
curve, typically by an order of magnitude. Nevertheless, 
the shape of the curves is similar to that given by the 
approximating curve. Singularity at T — > T c + demon- 
strated by the dashed line (see insertion in Fig. [4]) is due 
to the divergence of f a in this limiting case. 

The bulk viscosity ([T6|) includes a rather complicated 
factor Q a depending on the EOS used. Using the energy 
conservation for a system with temperature-independent 
masses of particles, one may present the result {jT6j) as 
follows: 



T 



d 3 p 
(2tt) 2 



P 

3E a 



C 2 s E a 



(21) 



For a single-component gas, this expression coincides ex- 
actly with the 25-years-old result of Gavin [38[ . 

Chakraborty and Kapusta (CK) [3^] presented another 



expression 



Cck 



45 



T 



d 3 p 



3 f a F?(l + F?)Qi (22) 



which differs from ([T6")) but also reduces to ([2~Tj) for 
m a = const. Note that they also disregard the QP inter- 
action term in the energy-momentum tensor; see Eq. (99) 
of their work. The reasons for the differences between 
(TIB")) and (12"2"[) are discussed in the Appendix. The dash- 
double-dotted line in Fig. @] demonstrates the Cchx/s 
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FIG. 4: The bulk viscosity to entropy density ratio for a 
glue matter. Solid and dashed lines are the results of our 
QP two-phase gluon-glueball model with two relaxation times 
as in Fig. [3] (the vertical short-dashed line joins boundary 
points of the mixed phase). The dash-double-dotted line is 
the calculation result with Eqs. p2[) and (|19[) . and the dot- 
dashed one is calculated according to Eqs. (|23p and (|20|l . 
The perturbative estimate (|24p is plotted by the dotted line. 
Experimental points are from [37| (empty squares) and [2^| 
(filled circles). The thin short-dashed curve corresponds to a 
simple approximating dependence (,/s = 0.02/ yjT/T c — 1 to 
guide the eye. 



ratio following Eq. (|22p with the relaxation time r g = f g 
given by Eq. flU]) (for f" 1 -> at T ->• T c + 0). We see 
that Eq. ([2"2")l yields a strong T suppression of the bulk 
viscosity at T Z 1.5T C , as compared to that given by Eq. 
(fTB)) (compare the dash-double-dotted and dashed lines 
in Fig. gj. 

A somewhat different expression for C w as used by 
Bluhm, Kampfer, and Redlich (BKR) jH). In their 
model, the bag constant B is a functional of the non- 
equilibrium distribution function. They found 



<{T)-T- 



dT 



(23) 



Here it was assumed that the QP interaction contributes 
to the energy-momentum tensor. Thereby, compared to 
(fTB)) there appeared the second term Tdm\(T)jdT in the 
square brackets of Eq. ([23)) . For constant masses, the lat- 
ter equation is also reduced to ([2"Tj) . Numerical calcula- 
tions with ([2"3")) (see the dot-dashed line in Fig. [T]) give rise 
to the Cbkr/s ratio, which dramatically decreases, being 
in large discrepancy with both the above models and the 
lattice data for T > 1.5T C . For T > 1.9T C , Eq. ^ 
becomes invalid, providing neg ative values. 
A perturbative estimate [37| gives 



(C/S)pert « 0.02^ 



(24) 



for 0.06 < a s < 0.3. Applying the T-dependent coupling 
constant for a s = g 2 (T)/A-K, we get a perturbative 
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estimate of the bulk viscosity to entropy density ratio 
(plotted by the dotted line in Fig. @|. As is seen, in 
the region of its applicability, Eq. (f24|) produces only 
slightly larger values of C/s than those given by our QP 
model. The new lattice GD calculations demonstrate a 
significant increase of the ratio C/ s at approaching the 
critical point (C/s ~ 0.5 -f- 2 at T = 1.02 T c ). These 
values are reproduced neither by our QP model using 
relaxation time Eq. (|19[) nor by the approximating short- 
dashed curve exploiting a simple T dependence of the C/s 
ratio. 

Comparing results presented in Figs. [3] and 21 we see 
that in the gluon phase in a narrow vicinity of the critical 
point [for (T - T c )/T c < 0.1], the ratio C/v > 0.1 reaches 
the value (/rj ~ 0.3 for T ->• T c + 0. The ratio sharply de- 
creases with an increase of the temperature up to values 
e/77 ~ 10~ 2 -7- 10~ 3 for T > 2T C . The smaUness of this 
ratio controls the violation of the conformal symmetry. 



V. CONCLUSIONS 

A quasiparticle approach has been applied to the 
SU (3) glue matter with temperature-dependent masses. 
Matching the pure gluon and glucball phase descriptions 
by means of the Gibbs conditions allows one to describe 
successfully this system in a thermodynamically consis- 
tent way both above and below the critical temperature 
T c . For thermodynamic characteristics, the quasiparti- 
cle model results are in good agreement with the latest 
lattice data. 

The constructed equation of state was used to calcu- 
late the shear and bulk viscosities in the relaxation-time 
approximation in a wide temperature range. The magni- 
tudes of the shear and bulk viscosities are mainly deter- 
mined by the value of the relaxation time. We exploited 
two different values for the relaxation time that have been 
used in the literature; see Eqs. ([TO]) and (|2"0|). With the 
relaxation time (|19p . the shear viscosity to entropy den- 
sity ratio rj/s fits rather well the scant lattice data. We 
found that the ratio rj/s undergoes a discontinuity at the 
critical temperature T = T c . At T slightly above T c the 
ratio rj/s has a minimum, the value of which is close to 
the AdS/CFT bound l/4n. Then 77/s increases with the 
subsequent rise of the temperature. The bulk viscosity to 
entropy density ratio C/s also has a break at T c . Then it 
monotonically decreases with the temperature increase. 
Although the calculated C/s ratio essentially underesti- 
mates the upper limits given by the corresponding lattice 
data, its temperature dependence is well described. 

Within our model, the ratio C/v — 0.3 at T — > T c + 
and it sharply decreases with the rising temperature up 
to values C/v ~ 10~ 2 + 10~ 3 for T > 2T C . 

We point out that although our QP model describes 
well the thermodynamical characteristics calculated on 
the lattice it does not guarantee that all important 
physics is incorporated, especially in a description of the 
vicinity of the critical point, where fluctuation effects are 



increased. Using the averaged value for the relaxation 
time can be considered only as a very rough approxima- 
tion. The dependence of the relaxation time on the cou- 
pling constant g is also not well defined. Simplifying, all 
existing approaches to the evaluation of the bulk viscos- 
ity use some Ansatze reductions yielding different results. 
We used the kinetic approach with the collisional source 
of the viscosity disregarding other possible sources [l(| ■ 
One such source is associated with the presence of a soft 
mode @ in the vicinity of the second and weak first-order 
phase-transition critical points. Also, since statistical er- 
ror bars are very large, new more certain lattice data are 
required to draw a more definite conclusion on the agree- 
ment or disagreement of the calculated C/s ratio with the 
lattice results. 

Acknowledgement s 

We are thankful to D. Blaschke, W. Cassing, M. Goren- 
stein, Yu. Ivanov, E. Kolomeitsev, N. Koshelev and 
V. Skokov for useful discussions. This work was 
supported in part by the DFG grant WA 431/8-1, 
RFFI grants 08-02-01003-a, Ukrainian-RFFI grant 09- 
02-90423-ukr-f-a, and a Heisenberg-Landau grant. 



VI. APPENDIX 

Deriving kinetic coefficients, the authors of j22l - [25| 
used the relaxation-time approximation to the kinetic 
equation presenting the collision integral as 

StF a = -8F a /r a [F loc -^}, 



(25) 



where 



5F a = F a (E a [F] ) - F^iE a [F l ° c -^}), 



(26) 

see, e.g., Eqs. (38) and (40) in (25[. Here it was 
assumed that the collision term should be zero for 
the global and local equilibrium states, i.e., for F a = 

F loc.eq.(j5 o [ F loc.eq.])_ 

Then after setting F a = F a occq - (E a [F loc ' cq -] ) on the 
left-hand side of the kinetic equation, one finds 



Tn[Jr loc.cq.] M dFl°c.cq.(ff a[F loc 
£r F loc.eq.]Pa q x » 



(27) 



see Eq. (2.3) in [22j and Eq. (42) in |25f. We stress 
that all quantities on the right-hand side of this equation 
including the relaxation time t are expressed in terms of 
the local equilibrium distribution functions. 

To derive expressions for the shear and bulk viscosities 
(TTB"]) and (|T6|) . one presents spatial components of the 
variation of the energy-momentum tensor of quasiparti- 
cles (JTDJ) as 



5T lk = 



<1T 



-^lSF a - V / dT^P^ 

E a [F] <W E 2 a [F io ^-} 



5E a 



E / dT 



PaPa 



E a [F io 



c.cq. 



;8F n . 



(28) 
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To avoid cumbersome expressions we omitted antipar- 
ticle terms. The reduction done in the second line in 
Eq. (|28[l is actually an Ansatz: we vary only the distri- 
bution function and do not vary quantities that depend 
on the distribution function implicitly (through phase- 
space integrals incorporating the distribution function), 
i.e., SE a [SF] are set zero. This reduction is in the spirit 
of the relaxation-time approximation to the kinetic equa- 
tion, where the momentum-dependent relaxation-time 
parameter is replaced in actual calculations by an aver- 
aged value. Note that dropping the SE term, we actually 
ignore a sub-leading term in the case of a weak coupling 
constant and/or for a very dilute system, see Eq. (5.22) 
of [10]. The distribution function SF a counted from the 
local equilibrium value enters the expression for 5T lk , the 
expression for i5T 00 = [see (45) and (50) in [25|], and ex- 
pressions for the viscosities. Thus one can easily compute 
kinetic coefficients knowing thermodynamic quantities in 
the local rest frame u = 0. 

However, we should note that in the QP Fermi liq- 
uid theory following the work of Abrikosov and Kha- 
latnikov [41|, one usually uses a different procedure to 
obtain transport coefficients; see [42j for details. One 
exploits that in the original Landau collision term, the 
following combination enters: 

5 lj2E a [F]) {F 1 (E 1 )F 2 (E 2 )(lTF 3 (E 3 ))(l T F i (E 4 )) 



-F 3 (E 3 )F 4 (E 4 )(1 T T F 2 (E 2 ))} , 



where E a are functionals of the exact non- equilibrium dis- 
tribution function, a = 1,2,3,4. The term in curly brack- 
ets is zero not only for F a = i^oc.eq.^^loc.eq.)] but 

also for F l ° c - cc i-(E a [F}). Thereby, StF l a oc - cc iiE a [F]) = 0. 
Thus introducing 



SF a = F a (E a [F}) - F l ° c ^(E a [F}) 



(30) 



in the relaxation-time approximation, we may rewrite the 
collision term as 



StF a = —5F a /T a [E(F)}. 



(31) 



The quantity f[F] entering Eqs. (f3"T|) and ([22"]) depends 
on the unknown exact non-equilibrium distribution func- 
tion, since the (5-function term in the collision integral 
and the local equilibrium distributions there continue to 
depend on exact energies in this approach. If we want to 
calculate the value of r[i^ loc oq ] entering Eqs. (JT5J) and 
(fT6|) using Eq. ([31]), we should still expand E[F) in (|3"T]) 
near the known value E[F loc ' cq '] everywhere including 
the 5-function term in the collision integral. 

From the left-hand side of the kinetic equation, one 
gets 



SF a 



f a [F] dFa° c - cq -(E a [F]) 



E a [F] dxZ 
Then one may use a simple expression for ST tk , 



ST 



E 



dT p^_ SF 



E„[F] 



(32) 



(33) 



since now variations are always performed at fixed E a . 
As above, the QP interaction term is omitted. Compar- 
ing the second line of Eqs. (|28|) and (|33|) . we see that 
disregarding the implicit dependence E[5F], Refs. [22| — 
[25l | actually do not distinguish between distributions SF 
and SF. 

Then in both considered approaches, one uses the exact 
relation E a [F] = ST 00 /SF a , i.e., that the variation of the 
energy is determined through SF a as 

ST oo = W dTE a [F]SF a 

a J 

= E / drE a [F loc -^]SF a + 0[(5F) 2 }. (34) 



Following (J26J) and <J30J), we have 



SF -SF = 



^^loc.cq.^^loc.cq.j-j 

dE 



SE. 



(35) 



Furthermore, instead of using a complicated implicit de- 
pendence SE[8F] with SF given by E q. |2"T1) . which would 
be the fully correct procedure, Ref. [39| uses the Ansatz 
relations [see Eq. (102) of that work] 



SF = exp 



^J^loc.cq.j ~l £;loc.cq 



-ST, 



(36) 



(29) 



SE 



ST mdm 
~~E dT 



loc.cq. 

SF T 2 mdm 
7£2 dT ' 



which assume that the distribution function in a nonequi- 
librium state has the form F = e~ E ^- F ^ T ^ in the Boltz- 
mann limit F « 1. Thus, although Ref. [39| distin- 
guishes between distributions SF and SF, it uses very 
special relations (|36[) . which might be incompatible with 
8E(m[SF}) as it follows from Eqs. ([32]) and ([TT]) . 

To find bulk viscosity, one further expresses F a = 
Fl oc - c i-(E a [F])(l - A a d p u p ), see Ref. [H], and one ob- 
serves that the shift of the solution A(E) — > A(E) — bE 
generates new solutions of the Landau kinetic equation 
for arbitrary constant b. Then one chooses b to explicitly 
fulfill the Landau-Lifshitz condition u^ST^ = 0. Note 
that this modification of the solution is quite unneces- 
sary provided one guarantees that the condition ST 00 = 
holds in the local rest frame. We have checked that this 
condition is satisfied in our QP model. 

Finally, within this approach one arrives at the expres- 
sion ([2"2"]) for the bulk viscosity, which is explicitly posi- 
tive definite, whereas positive definitcness of Eq. ([16]) is 
not seen explicitly. However, we stress once more that all 
quantities in H22\) still depend on exact energies, while the 
way in which the quantities in \lb}) depend on unknown 
exact distribution functions is hidden. Thus explicit pos- 
itive definiteness of expression ([2"2"]l for £ presents actu- 
ally only an apparent improvement. In any case, to use 
Eq. ([2"2"]) in practical calculations, where only equilib- 
rium quantities are known, one should replace E[F] by 
E[F loc - cq -]. 
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Moreover, we should stress that the values of the re- 
laxation time in (fT5)) and (f212j) are different. Since we do 
not perform complicated microscopic calculations of the 
relaxation time but only estimate its average value, we 
actually cannot determine whether (fl~6)) or (|22[) is more 
preferable, and we may use both of them. 

Note that Eq. ([25)1 is derived for a different model, 



where the QP interaction contributes to the energy- 
momentum tensor. Also, authors of [32j use different 
value for the relaxation time. 

Thus different Ansatze used in derivation of Eqs. (|16p . 
and (|^|) lead to different values of the bulk viscos- 
ity, as is shown in Fig. HI 
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